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Abstract. In this paper, we construct smooth forward Ricci flow evolutions 
of singular initial metrics resulting from rotationally symmetric neckpinches 
on § n + 1 , without performing an intervening surgery. In the restrictive context 
of rotational symmetry, this construction gives evidence in favor of Perelman's 
hope for a "canonically defined Ricci flow through singularities" . 
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1. Introduction 

Finite time singularity formation is in a sense a "generic" property of Ricci flow. 
For example, on any Ricmannian manifold where the maximum principle applies, a 
solution (M. n ,g(t)) whose scalar curvature at time t — is bounded from below by 
a positive constant r must become singular at or before the formal vanishing time 
Tform = n/2r. In some cases, e.g. if the curvature operator of the initial metric 
is sufficiently close to that of a round sphere, the entire manifold disappears in a 
global singularity. A beautiful example is the result of Brendle and Schoen that 
a compact manifold with pointwise 1 / 4-pinched sectional curvatures will shrink 
to a round point [5]. Under more general conditions, one expects formation of 
a local singularity at T < oo. Here, there exists an open set f2 of M such that 
limsup t y T |Rm(a;, t)\ < oo for points x € f2. 

Such behavior has long been strongly conjectured; see e.g. Hamilton's heuristic 
arguments [13, Section 3]. All known rigorous Riemannian (i.e. non-Kahlcr) exam- 
ples involve "necks" forming under certain symmetry hypotheses. Local singularity 
formation was first established by Simon on noncompact manifolds [26] . Neckpinch 
singularities for metrics on S™ +1 were studied by two of the authors [1, 2]. Gu 
and Zhu established existence of local Type-II (i.e. slowly forming) "degenerate 
neckpinch" singularities on S" +1 [12]. 

Continuing a solution of Ricci flow past a singular time T < oo has always 
been done by surgery. Hamilton introduced a surgery algorithm for compact 4- 
manifolds of positive isotropic curvature [15]. (Huisken and Sinestrari have a related 
surgery program for solutions of mean curvature flow [16].) Pcrclman developed 
a somewhat different surgery algorithm for compact 3-manifolds [22, 23, 24]. A 
solution of Ricci Flow with Surgery (rfs) is a sequence (M^, g(t) : < t < TjT) of 
maximal smooth solutions of Ricci flow such that at each discrete surgery time Tu = 
Tj~, v the smooth manifold (M.^, +1 , g(Tf7 +1 )^ is obtained from the singular limit 
(M^, <7(T^)) = (M^,lim 4/ , T + g(t)~) by a known topological-geometric modification. 
Each geometric modification depends on a number of choices (e.g. surgery scale, 
conformal factors and cut-off functions) that must be made carefully so that critical 
a priori estimates are preserved. The technical details of RFS for a 3-manifold are 
discussed extensively in the literature; see e.g. [18], [6], [21], and [30]. 

It is tempting to ask whether the choices made in surgery can somehow be elimi- 
nated. Indeed, Perelman conjectures [22, Section 13.2] that the following important 
and natural question has an affirmative answer: 

Question 1 (Perelman). Let gi(t) denote a smooth Ricci flow so- 
lution obtained by performing surgery on singular data (M 3 , g{T)) 
at a scale hi > 0. Does the sequence {M 3 ,gi(t) : T < t < T,} have 
a well-defined limit as hi \ ? 

He writes: "It is likely that by passing to the limit in this construction one would 
get a canonically defined Ricci flow through singularities, but at the moment I don't 
have a proof of that." If one regards a sequence {(M n , gi(T)) : i e N} of surgically 
modified initial data as a sequence of smooth approximations to irregular initial 
data (M. n ,g(T)), Question 1 may be regarded as a problem of showing that the 
Ricci flow system of PDE is well posed with respect to a particular regularization 
scheme. 
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Somewhat more generally, Perelman's wish for a "canonically defined Ricci flow 
through singularities" might be rephrased as follows: 

Question 2 (Perelman). Is it possible to flow directly out of a 
local singularity without arbitrary choices? 

In this paper, we consider SO(n+l, M)-invariant metrics on §™ +1 , and within this 
restricted context provide positive answers to both Questions 1 and 2 by exhibiting 
forward evolutions of the rotationally symmetric neckpinch. We say a smooth 
complete solution (M. n ,g(t) : T < t < T') of Ricci flow is a forward evolution of a 
singular metric (M n ,g(T)) if, as t \ T, the metric g(t)\o converges smoothly to 
g(T)\o on any open subset C M for which g is regular on 0. Thus we effectively 
let the Ricci flow pde perform surgery at scale zero. In so doing, we show that 
any forward evolution with the same symmetries as (M™ , g(T)) must have a precise 
asymptotic profile as it emerges from the singularity. Although the hypothesis of 
rotational symmetry is highly restrictive, 1 our construction provides examples of 
what a canonical flow through singularities would look like if Perelman's hope could 
be answered affirmatively in general. 

There exist a few examples in the literature of non-smooth initial data evolving 
by Ricci flow, none of which apply to the situation considered here. Berne lmans, 
Min-Oo, and Ruh applied Ricci flow to regularize C 2 initial metrics with bounded 
sectional curvatures [3]. Simon used Ricci flow modified by diffeomorphisms to 
evolve complete C° initial metrics that can be uniformly approximated by smooth 
metrics with bounded sectional curvatures [27]. Simon also evolved 3-dimensional 
metric spaces that arise as Gromov-Hausdorff limits of sequences of complete Rie- 
mannian manifolds of almost nonnegative curvatures whose diameters are bounded 
away from infinity and whose volumes are bounded away from zero [28]. Koch 
and Lamm demonstrated global existence and uniqueness for Ricci-DeTurck flow 
of initial data that are close to the Euclidean metric in L°°(R n ) [19]; this may be 
regarded as a generalization of a stability result of Schnurcr, Schulze, and Simon 
[25]. In recent work, Chen, Tian, and Zhang defined and studied weak solutions of 
Kahler-Ricci flow whose initial data are Kahler currents with bounded C 1,1 poten- 
tials [8]. The special case of conformal initial data e u g on a compact Riemannian 
surface (M 2 ,g) with e u e L°°(M 2 ) was considered earlier by Chen and Ding [7]. 
(Though their proofs are quite different, both of these papers take advantage of cir- 
cumstances in which Ricci flow reduces to a scalar evolution equation.) The term 
"weak solutions" was used in a different context by Bcssicrcs, Bcsson, Boileau, 
Maillot, and Porti to describe certain solutions of RFS on compact, irreducible 
non-spherical 3-manifolds [4] . In another direction, Topping analyzed Ricci flow of 
incomplete initial metrics on surfaces with Gaussian curvature bounded from above 
[31]. " 

This paper and its results are organized as follows. In Section 2, we use rota- 
tional symmetry to simplify our problem. The assumption of rotational symmetry 
generally allows one to reduce the full Ricci flow system to a scalar parabolic pde 
in one space dimension. For this problem, we were not able to find a convenient 
global description of the solution in terms of a one-dimensional scalar heat equa- 
tion. However, in Section 2, we show that, at least in appropriate local coordinates 



On the other hand, formal matched asymptotics for fully general neckpinches predict that 
every neckpinch is asymptotically rotationally symmetric [2, Section 3]. 
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valid in a neighborhood of the singular point, a forward evolution of Ricci flow 

<">. , £ = 

emerging from a rotationally symmetric neckpinch singularity at time zero is equiv- 
alent to a smooth positive solution of the quasilinear pde 

,~ „% i 9 n — I — v 2(n — 1) , 9N 

(1.2) V t = VV rr - \v r + «r + ^ 2 (v ~ V 2 ) 

emerging from singular initial data which satisfy 

(1.3) wwt(r) = [1 + o(l)]wo(r) as r \ 0, 
where 

(1-4) w<,(r) = 



\{n-l) 



— log r 

Note that a smooth forward evolution of (1.3) must satisfy lim t \^o v ( r , t) — v- ln it(r) 
at all points where the initial metric is nonsingular, i.e. at all r > 0. 

There are only two ways that the solution (1.1) can be complete. The first is 
that v satisfies the smooth boundary condition v(0, t) = 1 for all t > that it exists. 
Because lim. r ^o Winit( r ') — 0, this is incompatible with the initial data, meaning that 
v immediately jumps at the singular hypersurface {0} x S™, yielding a compact 
forward evolution that heals the singularity with a smooth n-ball. In this case, 
the sectional curvatures immediately become bounded in space, at least for a short 
time. The second possibility is that v remains singular at r = 0, but the distance to 
the singularity measured with respect to g(t) immediately becomes infinite, yielding 
a noncompact forward evolution, necessarily with unbounded sectional curvatures. 
In Section 3, we show that the second possibility cannot occur: we prove that 
any smooth complete rotationally symmetric forward Ricci flow evolution from 
a rotationally symmetric neckpinch is compact. (See Theorem 2.) This result 
contrasts with Topping's observation that flat IR 2 punctured at a single point can 
evolve under Ricci flow by immediately forming a noncompact smooth hyperbolic 
cusp, thereby pushing the "singularity" to infinity [31]. In this section, we also 
prove that any smooth forward evolution satisfies a unique asymptotic profile. (See 
Theorem 3.) 

In Section 4, we derive formal matched asymptotics for a solution emerging from 
a neckpinch singularity. This discussion is intended to motivate the rigorous argu- 
ments that follow. Because the initial lim,.\^ w init (r) = and boundary v(Q, t) = 1 
conditions are incompatible, one cannot have v(r,t) — > Vi n - lt (r) uniformly as t \ 0. 
The solution must resolve this incompatibility in layers. Consequently, we describe 
the asymptotic behavior of the formal solution for small t > by splitting the (r, t) 
plane into three regions: 



inner parabolic outer 



•log* 



Vt 



In Section 5, we make these asymptotics rigorous by constructing suitable sub- 
and super- solutions in each space-time region and ensuring that they overlap prop- 
erly. 
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Finally, in Section 6, we prove a compactness result that shows that a subse- 
quence of regularized solutions does converge to a smooth forward evolution from 
a neckpinch singularity. (See Lemmas 12 and 13.) 

Here is a slightly glossed summary of our results. (For more detail, including 
how the constants are chosen, see Lemmas 5 and 6 and Theorems 4, 5, and G.) 

Theorem 1. For n > 2, let (S n , <7o) denote a singular metric arising as the limit 
as t /" T of a rotationally symmetric neckpinch forming at time T = 0. Then 
there exists a complete smooth forward evolution (§ n+1 ,g(t) : Tq < t < T\) of go by 
Ricci flow. 

Any complete smooth forward evolution is compact and satisfies a unique asymp- 
totic profile as it emerges from the singularity. In a local coordinate < r < r, < 1 
such that the singularity occurs at r — and the metric is 



g(r,t) -- 

this asymptotic profile is as follows. 
Outer region: for c\\fi < r < C2, 



{drf 2 
v(r,t) +r 5can ' 



;(r,t) = [l + o(l)] 



n- 1 
—4 log r 



l + 2(n-l)- 



uniformly as t \ 0. 



Parabolic region: let p = rj\t and t = logi; then for C3/ 'v '— r < p < C4, one 



has 



v(r,t) = [l + o(l)] 



n — 1 

-2t 



2(n- 1) 



uniformly as t \ 0. 

Inner region: let a = \f—rp = W—r/t r; then for < a < 05, one has 
v(r,t) = [1 + o(l)]Q5( -) uniformly as t \ 0, 

where — r- + cr g can is me Bryant sohton metric. 

We recall the construction and some relevant properties of the Bryant soliton 
metric in Appendix C. 

Geometrically, the results above admit the following interpretation. In the inner 
layer, going backward in time, one sees a forward evolution emerging from cither 
side of a neckpinch cusp by forming a Bryant soliton, which is (up to homothety) 
the unique complete, rotationally symmetric steady gradient soliton on K n+1 . This 
behavior is unsurprising: as fixed points of Ricci flow modulo diffeomorphism and 
scaling, solitons are expected to provide natural models for its behavior near sin- 
gularities. Our asymptotics confirm this expectation and give precise information 
on the length and time scales on which the forward evolution is modeled by the 
Bryant soliton. 

It seems reasonable to expect that a solution will continue to exist if one admits 
initial data that are small, possibly asymmetric, perturbations of the data consid- 
ered here. It also seems reasonable that a uniqueness statement should hold, but 
proving that would likely require methods different from those used in this paper. 
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(Ml) 


ipo(s) > for all s S J 


(M2) 


Mo) = MQ = 


(M3) 


M*) = - 1 


(M4) 


V^) 2 = W + (i))r^ (-\o) 


(M5) 


^,( S )^ ( a ) = (2=l+ O (l))_i_ ( S \0) 


(M6) 


K(.s)|<l (0< S <£) 


(M7) 


3r # >o ^o( s ) ^ whenever ip (s) < 2r # 


(M8) 


3yi<ooV se j a (s)| < yi (where a (s) = ^qV'o ~ V'q 2 + *) 



Table 1. Assumptions on the initial metric go = (<is) 2 + ipo(s) 2 g t 



2. Recovering from a neckpinch singularity 

2.1. The initial metric and its regularizations. We will construct Ricci flow 
solutions starting from singular limit metrics resulting from rotationally symmetric 
neckpinches. For a review of such metrics, see Appendix A and, in particular, 
Lemma 14. 

Let NP and SP be the north and south poles of the sphere § n+1 . We identify the 
(doubly) punctured sphere §™ +1 \ {NP,SP} with (— 1, 1) x S™. On this punctured 
sphere, we then consider initial metrics g$ of the form 

(2.1) 50 = ^o(a;) 2 (^) 2 + V , o(^) 2 5can 

where fo,ipo ar e smooth functions on (—1,1). Table 1 lists the assumptions we 
make about the initial metric. (These assumptions are satisfied by the singular 
limits of the solutions studied in [2].) To make our assumptions geometric, we 
break gauge invariance by choosing distance s to the north pole as the preferred 
coordinate. In this coordinate, the initial metric then appears as 

(2.2) g = (ds) 2 + V>o(s) 2 ffca„ 

where £ is the distance between north and south poles in the metric go, and where 
V>o S C°°(J) n C 1 ( J), with J = (0, £). 

Since the initial metric go is singular at the north pole, the standard short-time 
existence theory for Ricci flow does not provide a solution. Because our initial 
metric is of the form (2.2) with V'o(s) ~ s\ log s| -1 / 2 for s — > 0, the volume of a ball 




Figure 1. The initial data 
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with radius s centered at the singular point is (C + o(l))s n+1 | logs) - ™/ 2 = o(s" +1 ). 
So no continuous metric g on S™ +1 exists for which 3c > : eg < g can < c~ x g. 
Therefore, neither the methods of Simon [27] nor of Koch-Lamm [19] can be used 
here to construct a solution. 

Instead, we will construct a solution by regularizing the metric in a small neigh- 
borhood of the north pole, yielding a smooth metric <? w for each small ui > 0, with 
9u ~~ * 9o as w \ 0- Each of these may be regarded as a "surgically modified" 
solution, obtained by replacing the singularity with a smooth (n + l)-ball. The 
short-time existence theorem then guarantees the existence of a solution g^it) on 
some time interval < t < T u , with g u as initial data. We will obtain a lower 
bound T > for T u and show that a subsequence of the g^ converges in C°° for 
all < t < T. 

The lower bound for T w follows easily from our previous characterization of 
singular solutions in [1]. Obtaining compactness of the solutions thus obtained 
turns out to be harder and will consume most of our efforts in Section 6. 

In the following description of the regularized initial metric g u , we will have to 
refer to a number of constants and functions which are properly defined in Sections 
4 and 5. We note that while we have phrased the description in terms of the 
geometrically defined quantities s, tp, ip s , it applies to metrics of the form (2.1). 

For given small tu > 0, we split the manifold S" +1 into two disjoint parts, one of 
which is the small neighborhood N w of the north pole in which tpo(s) < p*y/u). 

On § n+1 \ 7i u , we let our regularized metric g^ coincide with the original metric 
go. Within 3\f w , we let g^ be a metric of the form g u = (ds) 2 + ip u {s) g ceta , where 
ip u is a monotone function, at least when ip u (s) < p**Juj. In 3\f w , we can therefore 
choose r = V'w(s) as coordinate. In this coordinate, we require g^ to be of the form 

(dr) 2 

(2.3) g u = + r 2 5can , 
where v u (r) satisfies 

(2.4) v~(r, ui) < Vu(r) < v + (r, ui). 

Here w ± (-,w) are the sub- and super- solutions constructed in Section 5, and eval- 
uated at time t — uj. 

Recall that the scale-invariant difference of sectional curvatures, 

a = (L - K)i? (= ^ ss - 4> 2 S + 1) 

is uniformly bounded from above and below under Ricci flow [1, Corollary 3.2]. Af- 
ter possibly increasing the constant A, we may assume that our regularized metrics 
g w satisfy hypotheses (M1)-(M8) in Table 1, with the exception that near s = 
one has tpL( s ) ~ ¥ +1 instead of (M4) and (M5) 2 . Also, because we have modified 
the metric go near the north pole, the distance between north and south poles 
will not be exactly £, although we do have £ u = £ + oil) as ui \ 0. 
Solving Ricci flow starting from g u produces a family of metrics 

g u (t) = tp(x, t) 2 {dx) 2 + ip(x, t) 2 5can 
which evolve according to (A. 3). For these metrics, one will then have both 

1^.1 < 1 



2 This actually follows from the fact that tp'(s) = y/v(r,u) 
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and 

M = |^v„-^ + i| <a 

Therefore, again by (A. 3), 

(2.5) |(V> 2 ) t | = \2a + 2n(ip 2 s - 1)| < 2(A + n) 

is bounded. In particular, one has {ip 2 )t < 2(A + n) at critical points of x i— ► ip(x,t). 
Thus if one sets 



Jo — 



then condition (M7) implies that for < t < T , the function x i— > ip(x,t) has no 
critical points with ip(x,t) < r#. In particular, there can be no new neckpinch 
before t = T$. We have thus proved 

Lemma 1. T w > T for all uj > 0. 

Since s i— ► ip u> {s 1 1) is monotone when ^ < a function r i— > u(r, i) is defined 
for < r < r# for which one has 

(2.6) ^ a (a;,t)=u(r,t)=u(V>(a:,i),*) 

for all a; on the northern cap of the sphere where ip < r# holds. With respect to 
these local coordinates, one may write the metric (2.2) as 

(dr) 2 2 

.9 = —TY2 +T Scan- 

The sectional curvatures K and L defined in (A. 2) are computed in these coordi- 
nates by K = —uu r /r and L = (1 — u 2 )/r 2 . 

2.2. Evolution equations for u(r,t) and v(r, t) = u(r,t) 2 . To derive a differen- 
tial equation for u, we first recall that ip s satisfies 

V^ ss + (n- 1) - 2 

ip tp* 

On the other hand, from the definition (2.6) of it, we find by the chain rule that 

d t ip s = tp t u r + u t . 

Combining these two relations, we get it t = dtip s — ip r Ur in terms of ip and its s 
derivatives. To rewrite this in terms of u and its derivatives, we use the fact that 
when t is fixed, it follows from u = ip s that 

d_ _ d_ 
ds dr ' 

One then arrives at a parabolic PDE for u, to wit, 

, . / x 2 n — 1 — u 2 . .u(l — u 2 ) 

(2.7) Ut — u[uu r ) r — uu r H u r + (n — 1) ^ . 

The quantity v = u 2 satisfies a similar but slightly simpler equation, 

(2.8) v t = J[v], 
where 

(2.9) $[v] = vv rr - \v 2 r + n ~l~ V v r + 2i ^^v{\ - v). 



d t ip s = ip sss + (n-2)^ip ss + (n-l)~ 



u—, and hence ip ss — uu ri ipsss — u(uu r ) ? 
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In terms of v, we have 

(2.10) K=* L=I^. 

2r r z 

2.3. The shape of the singular initial data in the r variables. We now derive 
the asymptotic behavior of the unmodified function v{r, t) as r \ at time t — 0. 
First we estimate Ui n it(?~) = u(r, 0). By equation (A. 4) of Lemma 14 in Appendix A, 
the quantity r(s) = tp(s,0) satisfies 

(2.11) r=[\V^l + o{l)]^J= 

as the arc-length variable s \ 0. If this were an exact equation, one could solve it 
explicitly for s, yielding 



W 



n - 1 



V n — 1 \ 2r- 

where W denotes the Lambert-W (product-log) function, i.e. the inverse of y <— > ye y 
restricted to (0,oo). For our purposes, it suffices to notice that (2.11) implies that 
log r — [1 + o(l)] log s as s \ 0, hence that 

(2.12) a= . 2 + o(l) r^-logr 

— 1 J 

as r \ 0. Again by Lemma 14, it is permissible to differentiate (2.11), whereupon 
one finds in equation (A. 5) that the initial data u; n it(') — V'sOfO) satisfies 

(2-13) u init (r) = + o(1) ] 

2 yf- log r 

as r \ 0. Squaring, we conclude that v initially must satisfy 

n — 1 1 
(2.14) Winit ( r ) = (_ +o( i))__ (rN0 ). 

3. Possible complete smooth solutions 

In this paper, we construct smooth solutions g(t) of Ricci flow on the compact 
sphere §™ +1 . However, our initial data are only smooth on the punctured sphere 
§™ +1 \{NP}. Inspired by Topping's example in [31], one could also look for solutions 
on the punctured sphere which at all positive times t represent complete metrics. 
Such solutions must necessarily be unbounded. Below, we will show that such 
solutions do not exist for our initial metric. Consequently, all possible smooth 
solutions starting from our singular initial metric extend to the compact sphere. 
Their sectional curvatures K and L must be bounded; and in view of L = {1—v) /r 2 , 
we find that the quantity u(r, t) must satisfy 

(3.1) Hr,t)-l\ <svL P \L g{t) \r 2 . 
Thus 

(3.2) v(0,t) = l 

is the only relevant boundary condition at r = for our problem. 

Our existence proof of a solution starting from the singular initial metric involves 
the construction of a family of sub- and super- solutions v s , indexed by parameters 
e, S > 0. In this section, we show that the v function corresponding to any smooth 



4 r(logr) 2 ' 
A(n- 1) f 2 1 
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solution must lie between the sub- and super- solutions v e s to be constructed in 
Section 5. 

3.1. Lower barriers for v. We first show that any such solution remains positive 
for small r and t, with an estimate that proves the following result. 

Theorem 2. Any smooth complete rotationally symmetric forward Ricci flow evolu- 
tion from a rotationally symmetric neckpinch singularity is compact with v(0,t) = 1 
for all t > that it exists. 

Given A S (0, 1), define v\(r) for < r < 1 by v\(0) = and 

(3.3) flA(r)=|= *("-!) 1 ( <r<l). 

4 — log r 

For later use, we note that for r E (0, 1), one has 

X(n — 1) 1 

(3.4) (v x ) r - 

(3 - 5) (v\)rr- 4 ^ r 2(_ logr )3 r 2( logr )2 

Theorem 2 is essentially a corollary of the following observation, which shows 
that on a sufficiently short time interval, V\(r) — e(l + t) is a subsolution of (2.8) 
for all small e > 0. 

Lemma 2. Let v be any nonnegative smooth function satisfying 

(3.6) v t >?[v], v(r,0)>v init (r), 

where 3 is defined by (2.9), and v- ln i t satisfies (2.14). 

Then v >v\ on (0,ro] x [0,io] for some small positive ro and to, independent of 
the boundary condition v(0, •). 

We note that the following argument also allows v to be a piecewise smooth 
supersolution whose graph only has concave corners, such as the "glued" super- 
solutions we construct in Section 5. We will only write the proof for the smooth 
case. 

Proof. Since A < 1, we may fix < r 1 so small that v\(r) < min {vi n [ t (r), j} 
and (logr) 2 > (n — l)/4 for < r < ro. Then fix to > small enough so that 
v(r, t) exists for < t < to- By making to smaller if necessary, we may assume that 
v\(r ) < v(r ,t) for < i < to- 
Given e > 0, define 

(3.7) f e (r, t) = v(r, t) - v x {r) + e(l + t). 

We shall prove that f £ (r, t) > for < r < ro- The lemma will follow by letting e 
go to zero. 

To simplify the notation, we henceforth write f = f £ and v = v\. We also set 
H=X(n- l)/4. 

Observe that /(r, 0) > e for < r < ro and that both /(0, t) > e and /(ro, t) > e 
hold for < t < to- To obtain a contradiction, suppose that there exists a first 
time is (0, to] and a point f € (0,ro] such that f(f,t) = 0. Then by (3.6), one has 

(3.8) 0>/ t (f,t) = 3>(f, *)]+£. 

We claim that 5"[v(f,t)] > 0. This contradicts (3.8) and proves the lemma. 
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To prove the claim, observe that 

= f r {f,t) = v r (f,t) - v r (f) and < f rr (f, t) = v rr (f, t) - v rr (f, 1). 
Thus we obtain 

f 2 J [v(f, t}] = vf 2 v rr — - (fv r ) 2 + (n — 1 — v)fv r + 2(n — 1)(1 — v)v 
> vf 2 v rr — - (rv r ) 2 + (n — 1 — v)rv r + 2(n — 1)(1 — v)v 

V \ i o V 

n — 1 — 2v — 



(logr) 2 |_ 2 (log r) 2 

+ 2«{(n-l)(l-i;) + 



(-logr) 3 



2(logr) 2 1 (logr) 
> 0. 

Here we used the facts that n > 2 and < v < v < 1/4 at (f, i) . This proves the 
claim and hence the lemma. □ 

Now we can prove the main result of this subsection. 

Proof of Theorem 2. Let % ro denote the complement in §™ +1 of the neighborhood 
of the north pole NP e § n+1 in which -0in.it; < r o- For a smooth forward evolution, 
% ro is precompact for small positive time. So the only way a noncompact solution 
could develop is if a neighborhood of the singularity immediately became infinitely 
long. But for all r G (0, ro], Lemma 2 implies that 

1 2 



y/v(r, t) \/n - 1 
It follows that 

d g (t)(0,r ) - / ri —^ < oo, 
Jo \Jv{r,t) 

hence that the solution is compact. Because the sectional curvature of the metric 
g on planes tangent to {r} x S™ is L = (1 — v)/r 2 , the solution will be smooth only 
ifv(0,f) = l. □ 

3.2. Yet another maximum principle. In this section, we prove that the v 
function for any smooth, complete forward evolution with go as initial metric is 
trapped between the sub- and super- solutions which we will construct in Section 5. 
This implies the claims about the asymptotic profile of the solution in Theorem 1. 
We begin by establishing a suitable comparison principle. 

Lemma 3. Let v~ and v + be nonnegative sub- and super- solutions, respectively, 
of vt — $[v] ■ Assume that either v~ or v + satisfies 

(3.9) v rr <C, K=~ <C, L = < C 

2r r z 

for some constant C < oo on a compact space-time set S = [0,f] x [0, t\. 

If v~(f,t) < v + (f,t) holds for < t < i, and v~(r,0) < v + (r,0) holds for 
< r < f, then v~ < v + throughout S. 
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111 this lemma we assume that v are smooth sub- and super- solutions. However, 
the proof works without modifications in the case where v ± are piecewise smooth, 
where the graph of u~ only has convex corners, and the graph of v + only has concave 
corners. In our maximum principle arguments, we shall only evaluate v ± at "points 
of first contact with a given smooth solution" which are necessarily smooth points 
of v . Thus the hypothesis (3.9) is to be satisfied at all smooth points of v + or v~; 
and, in particular, we do not intend to interpret the second derivative v rr in (3.9) 
in the sense of distributions. 

Proof. We prove the Lemma assuming that v + satisfies (3.9). 
For A > to be chosen later and any a > 0, define 

(3.10) z = e- xt (v + -v~) + a. 

Then z > on the parabolic boundary of 3. We shall prove that z > in 3. 
Because this implies that v + — v~ > —ae xt in 3, the lemma follows by letting 
a\0. 

Suppose there exists a first time t £ (0,t| and a point r £ (0, f) such that 
z(r,t) = 0. Then z t (r,t) < 0, and 

v + (r, t) — v~(r, t) — ae xt , v^(r,t) — v~(r,t), v^ r (r, t) > v~ r (r, t). 

Hence at (r,t), one has > e xt z t , where 

e xt z t = w t + - V - X(v + - v~) 

( v + 1 

= v~[v+ r - v~ r ] + (v~ - v + ) I A - u+ + ^ - 2(n - 1) — 

> (v- — v + ) < A — H — - — 2(n - 1) = 

Thus using the uniform bounds v^ r < C, —v^/r < C, and (1 — v + )/r 2 < C on 3, 
together with w~ > 0, one obtains 

0>e xt z t >ae xt (X-C). 

This is a contradiction for any A > C. The result follows. 

To prove the Lemma in the case that the subsolution v~ satisfies (3.9) one uses 
the fact that at a first zero (r, t) of z one has 



e 



xt z t = v + [v+ - v~] + (v~ - v + ) { A - v~ + — - 2(n - 1) 



v~ , A 



r 



□ 



Note that a different formulation of the lemma above is: if v + is a supersolution 
of Vt = 3\v], then v + + ae xt is a strict supersolution. 

Theorem 3. Let v denote any solution of the Cauchy problem 

v t = 9[v], v(r,0) = v ini t(r), 

that is smooth for a short time < t < t\ . Let V e s denote the sub- and super- 
solutions, depending on s,6 > 0, that are constructed in Section 5. 
For all small enough s,S > 0, there exist f £ ^,i e ^ > such that 

(3.11) v- s (r,t)<v(r,t)<v+ s (r,t) 

for all (r,t) £ 5 = [0, f] x [0,t\. 
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Proof. Let s, 5 > be given. 

Since v assumes our initial data, we have v(r, 0) = (1 + o(l))vo(r) (r \ 0). Since 
the sub- and super- solutions are initially given by vf g (r, 0) = (1 ± S)vo(r), we find 
that there is some f > such that v~ 5 (r, 0) < v(r, 0) < vj s (r, 0) for all r £ (0, r]. 

The solution v and the sub- and super- solutions vf s are smooth for r > 0, so 
there is some t > for which v~ s (f, t) < v(r, t) < v^ s (f, t) holds for < t < t . After 
shrinking t if needed, we may also assume that (1 — S)vo(f) < v(r, i) < (1 + 5)vo(f) 
holds for < t < i. 

Lemma -i would now immediately provide us with the desired conclusion, but 
unfortunately neither v nor v^ s meet the requirements needed to apply that result. 
We overcome this problem by comparing time translates of v and vf s . 

First we show that v > v~ s . Translate the given solution in time by a small 
amount n > 0, i.e. consider the function 

v(r, t) = v(r, t + k) 

on the space-time domain S K = (0,r) x (0, i — k). By Lemma 2, we know that 
v(r,t) > v~ s (f,t) > (1 — S)vo(r) for < t < t . Since the time-translate v is smooth, 
it does satisfy (3.9), and we may conclude that v(r, t) > v~(r,t) on 3 K . Letting 
k \ 0, we end up with i>(r, t) > v~ (Y, t) on H. 

Next, we argue that v < v^ s on S. To this end we introduce 

v*(r,t) = v+ s (r,K + t). 

Since vf s is a supersolution that starts out with vf s (r, 0) = (1 + 5)v (r), Lemma 2 
tells us that v^ s (r,t) > (1 + 5)v (r) > v(r, 0) for all (r, t) £ S. We already have 
v^~ s (f, t) > (1 + S)v (r) > v(f, t) for < t < t. So, since v* is piecewise smooth and 
satisfies (3.9), we obtain v(r,t) < v*(r,t) on S K . Letting k \ then completes the 
proof by showing that v < v^ s on S. □ 

4. Formal matched asymptotics 

As a heuristic guide to what follows, we now construct an approximate solution 
of (2.8) that emerges from initial data satisfying v(r, 0) = V[ n it(r) — [1 + o(1)]vq(t) 
as r \ 0, where by (2.14), is given by (1.4). By Theorem 2, we may restrict 
our attention to smooth solutions satisfying the (incompatible) boundary condition 
v(0,t) = 1. 

To describe the asymptotic behavior of the formal solution for small t > 0, we 
split the (r, t) plane into three regions, which we label inner, parabolic, and outer, 
as in the introduction. 

We shall describe the solution in these separate but overlapping regions, working 
our way from the outer to the inner region. 

4.1. The outer region (r ~ 1). Away from r = 0, we expect the solution to be 
smooth. So a good approximation of the solution at small t > should be 



v(r,t) « v (r) + tm(r), 
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where v% = 2F[vq]- One computes that 



''o 



3\vo] = ~{2(n-l)(l-«o) + 



- log '/ 



(n - 1) - 2v a + 



3v 



—2 log r 



— 



2(n -l) + 



— log r 



So we set 

(4.1) v out (r,t) = v (r) 



1 + 2(n- 1)- 



n- 1 
-4 log r 



1 + 2(77-1)- 



This suggests new space and time variables 



(4.2) 



P = 



T = log t. 



With respect to p, which is bounded away from zero in the outer region, the outer 
approximation may be written in the form 



(4.3) 



v ut(r,t) = 



77- 1 

-2r 



1 



2(n - 1) 



+ 0(r- 2 ) 



as t \ 0, or, equivalently, as t \ — oo. 

4.2. The parabolic (intermediate) region (r ~ Vt). With p and r given by 
(4.2), we define W by 

W{p,-\ogt) W(p,r) 



(4.4) 



v(r, t) 



— log t — T 

Then it is straightforward to compute that W satisfies 



(4.5) 



1 



W T + — {W - Qpar[V^]} = H paT [W}, 



where XL par is the first-order linear operator 



(4.6) 



77—1 p 



and Q par is the quadratic form given by 



(4.7) 



Qpar [Ml = WW, 



1 



pp 2 p 



(W p f 



W n 



1 



2(77- 1) 



2(71-1) , 

P P 2 



As t \ one has r = log t — ► — oo. So if the limit lim t \^ W(p, r) exists, equation 
(4.5) leads one to expect it to be a function Wo(p) which satisfies £ pa r[Wo] = 0. 
To get a better approximate solution, we can add correction terms of the form 
Wi(p)/(— t) % and substitute in (4.5). In this way one finds a formal asymptotic 
expansion of the form 

Wi(p) , W 2 (p) , 



(4.8) 



W( Pi t) = W (p) + 



(-r) {- T y 



in which the Wj can be computed inductively from 



(4.9) 



^par[W ] = 0, H paI [W j+ l] = (j + l)Wj - J2 QpartWi, Wj-i]. 



i=0 



Here Q par [0, V'] i s t nc bilinear form corresponding to the quadratic form Q par defined 
above. 



MINIMALLY INVASIVE SURGERY FOR RICCI FLOW 



15 



We will not use this expansion beyond the lowest order term, but it did prompt 
us to look for the sub- and super- solutions which we find in Section 5.3. 

Taking only the lowest order term, our approximate solution in the parabolic 
region is 

W (p) 



(4.10) 



where Wo is a solution of £ pa r[Wo] = 0. The general solution of £ pa r[Wo] 
W (p) = c [l + 2(n — l)/p 2 ], which gives us 

2(n-l)" 



is 



v par 



(r,*) = 



Co 



1 



Matching with (4.3) for fixed p and t 
c = (n — l)/2. Thus we get 

n - 



(4.11) 
and 



W (p) = 



1 



tells us that the constant Co should be 
2(n-l)" 



1 



-'par 



(r,t) = 



l 



l 



2(n-l) 



-2t 

For small p, more precisely for p = 0(l/\/— t), one defines a new space variable 
cr = p\J —t in order to write the approximation 

(n — l) 2 n — 1 



(4.12) 



r(r,t) 



cr -2r 

At first glance, it may be surprising that the approximate solution in the in- 
termediate region is found by solving the first-order equation (4.6) rather than by 
finding a stationary solution of a parabolic equation. This is caused by the fact 
that the parabolic PDE i> t = ^[v] is degenerate when v = 0. 

If the solution v were approximately self-similar with parabolic scaling, then one 
would have v(y/ip, t) w U(p) for some self-similarly-expanding solution U. However, 
equation (4.3) shows that this is incompatible with the behavior of v near the "outer 
boundary" of the intermediate region. 

Nonetheless, we remark that self-similarly-expanding Ricci solitons U = VV 
do exist. These were first discovered by Bryant, in unpublished work. Each is a 
solution of the ode 



(4.13) 



U 2 U, 



pp 



n-l-U 2 
P 



1 



1-U Z )U = 



derived from (2.7), but each emerges from initial data corresponding to a singular 
conical metric, 

(dr) 2 



(4.14) 



9 



U 2 



+ r'g c 



where > 0. We prove this assertion in Appendix B. 

We shall see below that solutions to (4.13) with U(0) = 1 and U(oo) = do exist 
when the term \pU' in (4.13) is absent. These solutions correspond to the Bryant 
steady soliton. 

Existence of solutions emerging from (4.14) also follows from Simon's work [27], 
at least when is close to 1; but Simon does not require the hypothesis of 
SO(ra + 1, R) symmetry. 
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4.3. The inner (slowly changing) region (r - W(-logi)) . The formal 
solution in the parabolic scale r = 0(y/T — t) found above becomes singular as 
r — > 0; and, in particular, it does not satisfy the boundary condition at r = 0. 
Thus we look for a "boundary layer" at a smaller scale which will reconcile the 
incompatible initial v(r,0) = and boundary v(0,t) = 1 (t > Q) conditions. Our 
derivation of the formal solution in the parabolic region suggests that the smaller 
length scale should be y/t/(—logt). So we let 



(4.15) 




and define 

(4.16) v(r,t) = V(a,e). 
Then V satisfies the PDE 

(4.17) 99 t {9V e -aV (7 }=? in [V], 

where 3^ is obtained by replacing r-derivatives in 2f with er-derivatives, namely 

? in[V ] = vv aa - \{v a f + n ~\~ v v a + 2 -^^(v V 2 ). 

We will abuse notation and simply write J for 3^. 

We begin with the observation that 99 t = 0((— logt) -1 ) = o(l) for small t, so 
that the crudest approximation of (4.17) is simply the equation 3\V] — 0. This ODE 
admits a unique one-parameter family of complete solutions satisfying 1^(0) = and 
V(oo) = 0. These solutions are given by 

V q ((t) = Q3(fccr) (0 < k < oo), 

where 05 is the Bryant steady soliton, whose asymptotic behavior is 2}(er) = cr~ 2 + 
o(er~ 2 ) for a — > oo. These assertions are proved in Appendix C. 

Equation (4.17) suggests that this crude approximation is off by a term of order 
99 1 , which prompts us to look for approximate solutions of the form 

(4.18) V(a,t)=V (a) + 6e t Vi(a), 

Here Vo(cr) — %$(ka) contains an unspecified constant k whose value we will deter- 
mine later by matching with the approximate solution from the parabolic region. 

To find an equation for V\, we observe that the LHS and RHS of (4.17) applied 
to V yield 

(4.19) 99 t (9V e - aV a ) = -96 t aV£(a) + (99 t ) 2 [Vi{a) - aV{(a)} + 9*9 tt Vi{a), 
and 

(4.20) 3\Vq + 99 t V x \ = 3\V ] + d3 Vo [96 t Vx] + o{00 t V{) = 99 t d? V(> [V,] + o{99 t V 1 ), 
respectively. Here d3 is the first variation of the nonlinear operator 5F, defined by 



de 
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It is given by the ordinary differential operator 



2 



(4.21) d? Vo =V (o-)~ 



n-l-V ((x) 



d 

da 



a 

+ r w _w + 2(»^i) 

(7 G 

We note that 

|0 3 tt | + |^ t | 2 = °W as i\0. 
So by keeping only the most significant terms in the lhs and RHS, we find the 
following equation for V%, 

(4.22) {dfvJiVx} = -aV£(a). 

We will first solve this equation in the case fc — 1 (when Vq — < B). The general case 
then easily follows by rescaling. 

Lemma 4. The ordinary differential equation 

d3<8{£} = -<t®'(ct) 

has a strictly positive solution £ : (0, oo) — > i/iai satisfies 

Ma 2 + o(a 2 ) (a\0), 
n-l + o(l) (cr/oo), 

/or some constant M > 0. 

^4iZ other solutions of d^^lC] = — <72}'(o~) i/iai are bounded at a = are given by 
C(cr) = C(o-) - Act£'(ct) /or an arbitrary X E R. 

IfVo(o-) = Q3(fccr) /or any fc > 0, i/ien 

£(fco-) 



(4.23) €(a) 



Vi(<r) = 



fc 2 



is a solution of (4-22). 



Proof. We know that ^[^(/co")] = for every fc > 0; differentiating this equa- 
tion with respect to fc and setting fc = 1, we find that the homogeneous equation 
dSFtg [(f)] — has a solution 

0(a) = -a<B'(a). 

Given that </> satisfies d7<% [</>] = 0, one can use the method of reduction of order to 
find a second (linearly independent) solution <j). One finds the following asymptotic 
behavior of (f> for small and large a: 

fCo-("- 1 ) + o(o--(™- 1 )) (<t\0), 
exp{-o- 2 /2(n-l) + o(o- 2 )} (a / oo). 



<f>(a) 



If € p is any particular solution of ^^[Cp] = — ct*B'((7), then the general solution 
to d9^[C] = -ctQ5'(ct) is 

(4.24) € g (a) = a4>{a) + b4>{<T) + <t p {a). 

To obtain a particular solution which is bounded at a = 0, we note that for small 
cr, the equation [£] = — erQS' is to leading order 

£" + _ 2{ - n - l) £ « -S"(0)a 2 , 

(7 cr 2 
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where Q5"(0) < by Lemma 18 in Appendix C. From this, one finds that a solution 



C p exists for which 



(4.25) € p (a) = (K + o(l))a 4 as a \ 0, where K = - 1 ' > 0. 

Since <j> is not bounded at a — 0, the only solutions given in (4.24) which are 
bounded at a — are those for which 6 = 0. 

Near a — oo, the equation d5F<B[C] = — is, to leading order, 

1 2(n- IK 2 

One then finds that there also is a solution (t^, which satisfies 

£00(17) = n—l + o(l) (<T — > 00). 

But £ 00 (cr) — £ p ((t), being the difference of two particular solutions, is a linear 
combination of and <f>. Since 0(er) = o(l) and <^(cx) = o(l) as it — > cx^, we 
conclude that one also has 

(4.26) €p(a) = n - 1 + o(l) 

as a — > 00. So we see that the general solution of GfffslC] = — crQ3'(cr) which is 
bounded as er — > is given by £(ct) = £ p (cr) + a<f>(a). Setting a = —A leads to the 
general solution described in the statement of the lemma. 

The particular solution £ p (<t) is positive for small and for large a. The solution 
—<p(cr) = — cr!B'(cr) to the homogeneous equation is positive for all a > 0. Hence if 
one chooses A sufficiently large the resulting solution £ = £ p (cr) — \o-*B'(cr) will be 
strictly positive for all a > 0. This is the positive solution of dCF^C] = — a*B'(a) 
which was promised in the lemma. 

Finally, if d£F<8[(£] = — <j*B'(cr) then one verifies by direct substitution that 
V (a) = «B(fccr) and Vi(a) = k- 2 <L{ka) satisfy (4.22). □ 

We return our attention to the approximate solution V(a, t) in (4.18). Combining 
Vo(cr) = (1 + o(l))(ka)~ 2 and V\{a) = n— l + o(l) for large a with the observation 
that 



1 -logf _ | + o(l) 



(*\0), 



V(a, t) = V (a) + OOtV^a) = [1 + o(l)] 



2(logt)2 -log* 

one sees that for large a and small t our approximate inner solution satisfies 

1 1 

fc^2 + _2( n - 1) k 2 log t 

If we try to match this with the "small p expansion" (4.12) for our approximate 
solution in the parabolic region, then we see that we should choose 

(4.27) k = ' 



n — 1 

5. Construction of the barriers 

5.1. Outline of the construction. In this section, we will construct lower and 
upper barriers for the parabolic PDE 

v t = ^[v] = <vr 2 v rr — - (rv r ) 2 + (n — 1 — v)rv r + 2(n — 1)(1 — v)v \ . 
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Outer 


(l±6)v (r) + (l±e)t?{(l±S)v (r)} 


p*Vt < r < 


Parabolic 


— T T^p* 


J— < p < 3p* 


Inner 


<8(fc±cr) + (1 T s)66 t k 2 ± £(k ± a) 


< a < 3cr, 



Table 2. The sub- and super- solutions with their domains. Here, 
8 and e are sufficiently small; = Ae^ 1 / 2 , = Be^ 1 ^ 2 for certain 
constants A, B depending only on n; and j±, k± are given by (5.2) 
and (5.3). 



These barriers will apply to initial data satisfying 

Uinit(r) = [1 + o(l)]v (r) as r \ 0, 
where vq is the asymptotic approximation defined in (1.4), namely 

n- 1 



—4 log r 

The barriers will be valid on a sufficiently small space-time region 

< r < r*, 0<£<£*. 

Note that r* will not exceed the quantity r# from assumption (M7) concerning the 
initial metric. (See Section 2.1.) 

Because we will not be able to write down barriers that are defined on this whole 
domain, our construction proceeds in two steps. Theorems 4, 5, and 6 constitute 
the first step. In this step, in accordance with the matched asymptotic description 
of the solution in Section 4, we will produce three sets of barriers, each in its own 
domain. (See Table 2.) Note that the domains overlap. In all three cases, time is 
restricted to < t < i*. The parameters < r#, p*, er*, and will be defined 
during the construction. Although the construction admits free parameters 7, 6, e, 
and k, all but 8 and e will be fixed in the second ("gluing") step. 

After constructing separate barriers, we must "glue" them together in order to 
make one pair of sub-/super- solutions. For example, to glue the subsolutions in 
the parabolic and outer regions, we define 



v (r, t) = 

and 



t (r, t) 3p*\/i < r < 

t (r,t) a*\J-t/\ogt < r < p*Vi 



v (r,t) = ma,x{v out (r,t),v pm .(r,t)} 

in the overlap between the outer and parabolic regions, when p^^/t < r < 3p*%/t. To 
be sure that this construction yields a true subsolution, we will verify the following 
"gluing condition" : 

(5.1) Vpax(p*Vi,t) > v out (p*Vi,t) and v pal (3p*Vt,t) < u ou t(3p* >/*)*)• 
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Lemmas 5 and G constitute the second step of the construction. When this step is 
completed, we will have chosen 

(l-S)e 2 1 1 

(5-2) 7-=*+ , \"1 r , and fc_ - 



e + 2A 2 /(n-l) ' n-l y/\ - 1 _ - £ /2(n - l) 2 

for subsolutions, and 

(5.3) 1+ = S + , i!,tf, )£2 T\ and k +- ' ' 



e + 2Ay(n-l) ^ n-l ^1 + 7+ - e/2(n - l) 2 

for supersolutions. Note that for small S, e > 0, one has fc_ > k + . Because 2J(cr) is 
decreasing, this implies in particular that V|~ < holds in the inner region. 3 

5.2. Barriers in the outer region. In Section 4.1, we constructed an approx- 
imate solution of the form v(r,t) — vo(r) + tui(r). It turns out that a slight 
modification of this approximate solution yields both sub- and super- solutions in 
the outer region. 

Theorem 4 (Outer Region). There exist positive constants A = A(n), = r*(n), 
and £* = f*(ro, e) such that for all \S\ < \ and for all small e > 0, the functions 

i4 t (r, t) = (1 ± 5)v (r) + (1 ± e)t?[(l ± S)v (r)] 

are sub- (v~ ut ) and super- (v^ ut ) solutions in the region 

fiout =§={(r,t):p,Vt<r<r„ Q <<<£*}, 

where — A/y/e. 

The constant (5 need not be positive in this theorem, however, if one wants a 
properly ordered pair of sub- and super- solutions, i.e. if one wants v~ nt < v^ ut , 
then one must choose S > 0. 



3 By taking e > sufficiently small, depending on 8 > 0, one can make > (n — 1) 1 > fc+; 
but we will not need this fact. 
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Proof. We will show that v~ ut is a subsolution in a domain p^\ft < r < r„, where 
r* < r# and » 1 are to be chosen. The proof that v^ ut is a supersolution is 
entirely analogous. 

To simplify notation, we define constants 

a =(1-5)2=1, b=l-e, 

and functions 

A(r) = — — , r(r) = ?[aA(r)]. 
— log r 

We henceforth write v~ ut (r, i) as 

(5.4) Ua , b (r,t)#aA(r) + 6tr(r). 

Observe that 

T(r) = ^ {2(n - 1) + (n - 1)(1 - 2a)A - 2aA 2 + |aA 3 } . 
It follows that there exists > such that 

r(r)>C^ for 0<r<r„. 

Because d t (v a ,b) = bT, we have 

d t {v a ,b) - SR.fc] = (b ~ l)r - {?[aA + far] - J[aA]} . 

Here the first term (6 — 1)T has a sign for < r < r*. The other term vanishes for 
f = 0. So for small t, continuity of dt{v a ,b) — ^[v^b] implies that dt(v a ,b) — 3\v ai b] < 
if b < 1. However, the size of the time interval on which this holds will depend 
on r and, in particular, will shrink to zero as r — > 0. Wc now make this argument 
quantitative. 

Using the splitting of 5F into linear and quadratic parts, one finds that 

9\aA + btT] - 9\aA] = btL[T] + Q(aA + btT) - Q(aA) 

= bt{L[T] + Q(2aA + btT, T)} 

where the bilinear form Q(-, •) is defined by polarization from the quadratic form 
Q(-), i.e. via 

Q(f,9) = l{Q(f + g)-Q(f-g)}- 

Both L[f] and Q(f,g) are conveniently estimated in terms of the pointwise semi- 
norm 

[fh(r) = \f(r)\+r\f r (r)\+r 2 \f rr (r)\. 
Indeed, for < r < r*, 

|£/l<^[/] 2 and \Q(f,g)\<^[f} 2 [9h. 
So for all r < r*, one has 

[aA] 2 <CA and [T} 2 < C^. 

Therefore on the region p*\fi < r < r* , one has 

*[!•], <|A. 
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For > 1, this implies that 

\L\T] + Q(2aA + Wr, T)| < CAr~ 4 + CA 2 r~ 4 + CA 2 r - 4 p- 2 

< CAr~ 4 . 

Observing that T > CAr~ 2 for < r < r*, we thus estimate 

t\L[T} + Q(2aA + btT, T)\ < %^ < %T. 

pi • pi 

Returning to our estimate for dt(v a ,b) ~ ^[va^b], we now have 

dt(v a , b ) - J[v a , b ] < (6_i)r+-^r. 

pi 

Because T > 0, we may conclude that v a ^ is a subsolution in the outer region 
p* yt < r < provided that 

p* >C{l-b)- 1/2 . 

Because 1 — b = e > 0, the theorem follows by taking A= C. □ 

5.3. Barriers in the parabolic region. In the parabolic region, we use the sim- 
ilarity variables p, t, and W defined in (4.2) and (4.4). According to (4.5) the 
function W satisfies D pa r[W] = 0, where 

£par[Wl = W T - £ par [V^] + —{W - Q par [^]}. 

— T 

Here £ par is as in (4.G) while Q par [H^] is the same quadratic differential polynomial 
as Q[v], but with all r-derivatives replaced by p-derivatives, namely 
1 

^2 



Q par [W] = ^{Wp 2 W pp -\(pW p ) 2 - pW p W -2(n-l)W 2 }. 



p* 

Theorem 5 (Parabolic region). Let e > 0, p* = A/y/e^ be as in Theorem 4- 

There exist B — B{n 1 e) > 1, and t* — £*(n, e) such that for any 7 with < 7 < 
\ , the functions 

are sub- (v~ ar ) and super- (v^ ar ) solutions in the region 



^par 



{(r,t):^ s] [^<r<3 P *Vt, < t < u} • 



Here ^ 

W±AP,r) = (1 ± 7 )Wo(p) ± -J^, 
where Wo(p) is as in (4. 11). 

Proof. We consider the case of a subsolution. Recall that cWo is the general solution 
of the first-order ode £ par [VK] = 0. Moreover, 

(5.5) £ par [p- 4 ] = -2p" 6 (n - 1 + p 2 ) < -2(n - l)p~ 6 . 

To show that W~ ar is a subsolution, we will verify that DparfWj^J < on f2 par . To 
simplify notation, we write W = W~ ar for the remainder of the proof. 
Observe that 

[W] 2 = \W\ + p\W p \ + p 2 \W pp \ < C (p- 2 + S 2 |r|-V- 4 ) 



MINIMALLY INVASIVE SURGERY FOR RICCI FLOW 



2.i 



if p < 3p* = 3A/^fe, where C = C(n,A,e). In the parabolic region f2 par defined 
above, one has 

(5.6) < -. [ 2 < e, 



whence we get 
and also 

Now we compute that 



[W] 2 < Cp- 



(-T)V paI [W] = -J—r^ + B2 ^{P~ 4 } + (1 - 7)^0 - Qpar[W]. 



2B 2 

Assuming that r < —1, we find, using (5.5), (5.6), and also p < 3p* = 3Aj ^fe, that 

C-2(n- l)B 2 



(-r)D par [W/] < 



P 6 



in the parabolic region. Hence we conclude that the function W = W~ av will indeed 
be a subsolution provided that B 2 > max{l, C/2(n — 1)}. 

The left- and right- end points of the parabolic region at any time r are given 
by p = B I \l s(—t) and p = 3p* = 3A/^/e, respectively. So this region will be 
nonempty if — oo < r < -B 2 /(3A) 2 . Thus we choose = - max { 1, B 2 /(3A) 2 }. 

Construction of supersolutions W^ al . is similar. □ 

5.4. Barriers in the inner region. In the inner region, we work with the space 
and time variables a, 8 defined in (4.15). We consider V(a, 9) = v(r, t), as in (4.16). 
Then, according to (4.17), Ricci flow is equivalent to T>i n [V] = 0, where 

T>in[V] =06 t (6Vg - aV a ) - ?[V] 

= 6 2 V t - 69 t aV a - J[V]. 

The formal solution we found in Section 4.3 is of the form 

V(a,t) = V (a) + X69 t V 1 (a), 

where 

V (a) = S(fco-), and Vi(a) = fc- 2 C(fccr). 

In Section 4.3, we chose A = 1 and k — l/(n — 1) in order to match this solution 
with the formal solution in the parabolic region. Here we will show that small 
variations in k and A lead to sub- and super- solutions. 

Theorem 6 (Inner Region). Let e and B be as before. There exists t* = t*(n, e, B) 
such that for any k G [ 2(^-1) ) n^i\> ^ e functions 

<M) = V±(a, 6) = <B(ka) + (1 T e)96 t k- 2 £(ka) 
are sub- (f;~) and super- (Ujjj) solutions in the region 



<> m = {(r,t) :0<r<3-|W-— , ■ / / 
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We draw the reader's attention to the fact that 99 t €(ko-) > 0. So for fixed k, 
the subsolution V~ is larger than the formal solution (which has e — 0), while the 
supersolution is smaller. To get a properly ordered pair of sub- and super- 
solutions, we must (and can) choose V^I and V-~ with different values of k. 

Proof. We will prove that 

Vr = V (a) + (1 + e)ee t Vi(a), with V (a) = 58(fea) and Vi(cr) = k 2 €(ka), 

is a subsolution in the region f2; n , i.e. for 

< cr < 3cr* = 3.B/Ve, < < < i*, 

where i* is suitably chosen. The proof that Vj^ - is a supersolution is similar. 
Upon substitution, we find that 

(5.7) D in [^-] = 2 (l+e)(^ t ) t y 1 -^ f( ry o '-(l+e)(^ t ) 2 ( Ty i '-?[Uo + (l+e)^F 1 ]. 

We can expand the last term, keeping in mind that ^[Vq] — 0, and that l J\V] is a 
quadratic polynomial in V and its derivatives. We get 

?[V + (1 + e)90tVi\ - (1 + e)^ t dJvb^i] + (1 + s) 2 (99 t ) 2 Q[V 1 ], 

where 

q[v] = - i(K) 2 - - ^iy 2 

is the quadratic part of 3\V]. Applying this expansion to (5.7), we find that 

(5.8) V in [Vr] = £99 t aV^ + {l+e)9 2 {99 t ) t V 1 -{l+e){99 t ) 2 aV[-{l+e) 2 {99 t ) 2 Q[V 1 ] 

The key to our argument is that the first term dominates the others on the interval 
< a < cr* = 3B/y/e. 

The facts from Lemma 18 that the Bryant soliton 2}(-) is strictly decreasing, and 
that it is given by 25(cr) = 1 + 62c 2 + • • • for small cr, with 62 < 0, tell us that there 
is a constant 77 > such that 

(5.9) - erVoV) > r]a 2 for a e (0, 3a*). 

The asymptotics of £(cr) both at a — and a = 00 from Lemma 4 tell us that for 
some C = C(n, £) < 00, one has 

(5.10) |Vi|+ \vV{\ + \Q[Vi]\ < Co 2 for all a £ (0,3(7*), 

provided that cr* > 1. Finally, by direct computation, one finds that 

11 H 

9dt = ^^bgl + FToW)' 92{e9t)t = 2(-logt)3' 

so that |6> 2 (6>6» t ) t | < C(99 t ) 3 for some C < 00, and for small t. 
Together with (5.9) and (5.10), we find that 

(5.11) D[F in ] < {-ve + C99 t } 99 t a 2 

for all a 6 (0, 3cr*). Since 99 t = o(l) as t \ 0, we find that V- m is indeed a 
supersolution for small enough t. □ 
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5.5. Gluing the outer and parabolic barriers. The barriers Wp aI constructed 
in Section 5.3 generate sub- and super- solutions 

W ± 

4- par 



par -log* 

for the original equation vt = ^[v] in the parabolic region. 
Lemma 5. Let A, e and 6 be as before, and set 

(1 - S)e 2 

( 5 - 12 ) 7 = e) = 5 + \ -. 

e + 2 A 1 1 [n — 1) 

If — r* is sufficiently large, then the gluing condition (5.1) is satisfied for all r < r*. 
Proof. We will verify the relations 

V pax (p*Vt, t) > Uout(P*VM) and Upai(3p*Vt,t) < Uout(3/3*VM) 

in the W(p,r) rather than the v(r,t) notation. 

When written in the (p, r) variables, the subsolutions from the outer region take 
the form 

where 

If P* < P < 3p», then 

A = — — = = -2T- 1 + 0(r~ 2 ). 

- log r - log p — | log t 

Hence 

W mt (P, r) = (1 - 5)^ [l + ^-y-^) } + 0(T -i), (r - -oo), 

uniformly in p* < p < 3/9*. 

The subsolutions from the parabolic region satisfy 

W~ ai (p,r) = (1 -7)^ {l + ^r 11 } + 0(0, (r - -oo). 
The outer and parabolic subsolutions have limits as r — > — oo, namely 



out vp; ■ v ^ 2 \ v P 2 

w'-(p)^d-i)^- 1 ^ 1 2(71 ~ ir 



P 2 

respectively. At p = 0,oo, one finds that 

(1-8 

Txr-oo/ \ as p — > oo , 

ft 1q s gout P J 1-7 

Wpar (p) aSp^O. 

^ 1 — 7 
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Requiring W ~?°(2p*) = W p -?°(2p*) leads to (5.12). (Use = A/y/e.) If (5.12) 
holds, then p — 2p* is the only solution of W~^(p) = W~^(p), and it follows from 

(5.13) that 

W m T(P*) < W-™(p*) and W- U t(3p*) > W£>(3p.). 

In particular, the gluing condition (5.1) is met. 

Because W~ ut and W~ w are small (0(r -1 )) perturbations of and W~£° , 

respectively, these inequalities will continue to hold for all sufficiently large — r. □ 

A similar statement holds true for supersolutions. 

5.6. Gluing the inner and parabolic barriers. Recall that the inner and par- 
abolic regions are 

in = {0 < a < 3a-*, 0<t<t*}, 
n par = {a* < a < 3AV=r/Ve, < t < u} , 
respectively, where 

o~ * = B/^fe. 

We now verify the "gluing condition" between the inner and parabolic regions. 
Lemma 6. // 

(5.14) k = k-(5,e) ' ' 



n-1 y/\ _ 7 - e /2(n - l) 2 
ofid z/i? is sufficiently large (depending only on n), then the gluing conditions 

(5.15) Vr ( a „T) > V^a^r) and V£(3(t„t) < V^a^r) 
are satisfied for all t < t*, provided that — t„ is sufficiently large. 

Proof. By Theorem G, the subsolutions in the inner region f2; n have a limit as t \ 0, 
equivalently, as r \ — oo. The limit is 

(5.16) Vr°°(a) = HmV^^t) = B(fer). 

Our asymptotic expansion of the Bryant soliton (Lemma 18) implies that 

(5.17) <B(/ca) = (ka)- 2 + ((fccr)" 4 ) , (cr -► oo). 

By Theorem 5, the subsolutions in the parabolic region f2 par also have limits at 
t = 0, t — — oo, namely 

(l-7)(n-l) 2 B 2 

(5.18) V» = limVa>M) = " ^ L - ^ 

Equations (5.16) and (5.18) show that to establish (5.15), it will suffice to prove 

(5.19) Vr°° (a* ) > V~<? (a* ) and V-" 00 (3a, ) < V~<? (3a, ) 

and then to choose — t» sufficiently large such that VJ~(a, , r) > V^j.(a*,r) and 
^~(3a*,r) < Vp~ r (3a»,T) are preserved for all r < r*. 

Without loss of generality, we may assume that ~(n — 1) _1 < fc < 2(n — 
Then by (5.17), there is a constant C„ < oo such that 

(5.20) \^ 2 y- m °° ~ k " 2 \ < C n°~ 2 
for all a > 1 and all k under consideration. 
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Hence we have 

^ {y-^r - y°) = ^ 2 - a - 7)(n - 1) 2 + (s 2 + ^cy- 2 , 

where < 1. In order to verify (5.19), we must find the sign of the lhs for 
a = ma*, with m = 1 or m — 3. Using a* = £?£~ 1//2 , we find that when a = ma*, 

^n°° ~ V p -?) = fc- 2 - (1 - j)(n - l) 2 + ^(1 + VC n B- 2 ). 

In particular, (5.19) will hold if this quantity is positive for m = 1 and negative 
for m — 3. We can achieve this by first choosing B so large that we can ignore the 
term containing d: B > B n = \/100C n will do. Then we choose k so as to satisfy 

fc- 2 -(l- 7 )(n-l) 2 = -|. 

Solving this for k leads to the value fc_(<5, e) mentioned in the Lemma. Once /c is 
given this value and B is chosen large enough, (5.19) will hold. 

□ 

A similar statement holds for supcrsolutions. 

6. SUBSEQUENTIAL CONVERGENCE OF THE REGULARIZED SOLUTIONS 

In this section, we prove compactness of the family {gui{t) | < t < T} for some 
< T < To. We find a convergent subsequence g Uj {t) — > g*(t); we show that g*(t) 
is a smooth solution on S" +1 for small t > 0; and we verify that g*(t) is indeed a 
forward evolution from the singular initial metric go. 

We do this directly, rather than by invoking Hamilton's compactness theorem, 
which instead gives ^ (5^ ) (*) — * ^ or some sequence : S™ +1 — * S" +1 } of 

time-independent diffeomorphisms fixing the north pole NP £ § n+1 [14]. 

The main problem in establishing compactness is that, although we have precise 
control of the v function near the singular point, this information is only valid in a 
neighborhood of the form {P £ § n+1 | ip(P, t) < f} for some f > 0. 

6.1. Splitting S" +1 into regular and singular parts. The following lemma 
allows us to split the manifold S" +1 into a regular and a singular part. 

Lemma 7 (Existence of collars). For any < m < 1 and small enough a > 0, 
there exists T = T(a) € (0,Tq) such that for < t < T{a), the function 



r — r 



2 ' 



v co \ + {r, t) = min |l,me 3t/a2 + 
is a supersolution of Vt = 3^[v], while 

v co i~(r, t) — max |o, me~ l t a — 
is a subsolution. 

Proof. The statements are verified by direct computation. For G±{z) — m±z 2 one 
computes that 

? \G+ I ] I = a- 2 [2m + 0(a)] , 
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and 



2f 




f r — r\ 


= a- 2 


2m 













2{ r —) 2 + 0{a) 



hold when |r — r| < a, and when a < |f. 

If we now let m vary with time, we find that m(t) — z 2 is a subsolution for small a 
if m'(t) < cm(t) /a 2 for any constant c < 0: we choose c = — 1 to get the subsolution 
in the Lemma. Similarly, m{t) + z 2 is a supersolution if m'(t) > cm(t)/a 2 for any 
constant c > 2; we choose c = 3. □ 

In Theorem 3, we saw that a smooth forward evolution (if one exists) must lie 
between the sub- and super- solutions constructed in Section 5. Now we show, using 
Lemma 7, that the regularized solutions g u (t) obey a similar bound, uniformly for 
small uj > 0. 

Lemma 8. Let e,S > be given. Then there exist f, i > such that for each 
solution gu>{t), with oj > sufficiently small, one has 

(6.1) v- s (r,t + u) <v u (r,t) <v+ s (r,t + cj) 
for all re (0, f) and t £ (0, i) . 

Proof. To apply our comparison principle, Lemma 3, we need to show that the 
given solutions v^(r, t) cannot cross the barriers v e s (r, f) at the right edge r = f of 
their domain. The barriers from Lemma 7 allow us to do this. 

We let f be as in the construction of the sub- and super- solutions vf s . Then we 
choose m± and a > so that the sub- and super- solution from Lemma 7 satisfy 

w coi( m -> a ; r >°) < u mit(r) < v+ ol (m +) a; r,Q) 
for all r, and also so that 

v~ s (r,0) < v~ oi (m-,a;r, 0) and v+ s (f, 0) > v+ ol (m + ,a;r,0). 
Finally, we choose t > so small that this last inequality persists for < t < i . □ 
6.2. Uniform curvature bounds for <? w (i) for t > 0. 

Lemma 9. The curvature of the regularized solutions g^i) of Ricci flow constructed 
in Section 2.1 is bounded uniformly by 

(6.2) |Rm| < C— ^ 

for < t < t*, where the constant C depends only on the initial metric go and 
chosen sub- and super- solutions v . 

Proof. For some small enough e,5 > 0, we consider the sub- and super- solutions 
v^r s constructed in Section 5, which are defined in the "singular region" tp < rv 

Here r% = r*(e,<5) <C r#. We now choose m* 1 and a < \r* so small that 
«c + d(^ *) = m+e 3t / a2 + (r - r*) 2 /a 2 > v init (r) 



and 



v coiM) = m e t/a - (r-r*) 2 /a 2 < v init (r). 



In the collar-shaped region \if) — r*| < \a 1 the quantity v is now a bounded solution 
of v t — ^[v] which is bounded away from zero. Interior estimates for quasilincar 
parabolic pde imply that all higher derivatives will be bounded in an narrower 
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collar \ip — | < | a. In the "regular region" which begins with ip > and extends 
to the south pole SP G S n+ , we can now apply the maximum principle to 

(6.3) <A|Rm| 2 + Cn |Rm| 3 

to conclude that |Rm| remains bounded for a short time there. Since all higher 
derivatives of Rm can be expressed in terms of v and its derivatives, we find that 
these too are bounded in the collar, and by inductively applying the maximum 
principle to equations for dt\ V f Rm| 2 — A|V £ Rm| 2 , one finds that higher derivatives 
of Rm also remain bounded in the regular region. 

It remains to prove that (6.2) holds in the singular region, namely the region 
where ip < r* . We do this in the next Lemma. □ 

Lemma 10. Let < v~{r, t) < v + (r,t) be sub- and super- solutions of x>t — J[v] 
which we have constructed on the domain < r < r*, < t < t*, and suppose 
that v is a solution of this equation which satisfies v~ < v < v + . Then there is a 
constant C such that 
1 - 



< C — for < r < r% and < t < U 
~ t 



The constant C depends only on the sub- and super- solutions v and the constant 
A from (M8) in Table 1. 

This Lemma implies that the sectional curvatures L = (1 — v)/r 2 and K = 
—v r /2r are uniformly bounded by C(— log i)/t for all solutions caught between our 
sub- and super- solutions. This completes the proof of Lemma 9 above. 

Note that the Lemma predicts that sup |Rc| ~ C(— log i)/i, which is not inte- 
grable in time. This is consistent with the fact that our initial metric go is not 
comparable to the standard metric on S n+1 . To wit, there is no constant c > 
such that cg S "+i < go < |5s»+ 1 - 

Proof. The supersolution satisfies v + < 1, so we immediately get v < 1. The 
subsolution is bounded from below by 

for some constant C. This implies that 

1 -v < 1 ~v~ < c ~log^ 

r 2 — r 1 — f 

To estimate the second derivative, we choose a small constant v > and split the 
rectangle £1 = (0,r„) x (0, t*) into two pieces, 

fti = \ (r, t) e : r 2 < v — * — I and Q, 2 = \ (r,t) € : r 2 > v — t - — 
I -log* J { -log* 

Throughout the whole region fi, the quantity a = r 2 (K — L) is bounded by (M8), 
Table 1, so that we get 

I =\K\<\L\ + \a,r 2 \<C^ + A 2 . 
2r t r z 

In the region f2 2) the second term is dominated by the first, and we therefore find 

that \v r /r\ < C(- log i)/t on tt 2 . 
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In the first region f2i, the equation v t — 3[v] is uniformly parabolic, so that 
standard interior estimates for parabolic equations [20, §V.3, Theorem 3.1] provide 
an (r, independent bound for v r . The following scaling argument shows that this 
leads to the stated upper bound for \v r /r\. 

Let (f, t) € r2i be given, and consider the rescaled function 

1 — v(f + fx, t + f 2 y) 
w(x,y) = . 

Then w satisfies 

, s dw , N , , 2(n - 1) , 

6.4 — = 1 - (3w)w xx + \l3w 2 x + , , w x + , M l - I3w)w. 
ay 1 + x (1 + x) A 

In the rectangle — \ < x < |, — \ < y <0, we have 

r = f(l + x) S (|f, |f) , 

and 

t > t = t + f 2 y > i- if 2 > f 1 - -g-) i > \t 



•log* 



if ** is small enough. 

Hence w will be bounded by 



C-log*_ 2 
< w < — — = — r . 



We now choose /? = C ~ p f 2 , which always satisfies P < Cv, since (f, *) € fit. 

Our function mi is a solution of (G.4) in which 1 — j3w = v > v~ is uniformly 
bounded from below, 1 + x > | is bounded from below, and for which we have 
found < w < 1. Standard interior estimates now imply that |i/j x (0, 0)| < C for 
some universal constant. After scaling back, we then find that 



Mr,*) | 

as claimed. 



~«>x(0,0) 
r 



-log* _ 



□ 



6.3. Uniform curvature bounds at all times, away from the singularity. 

The initial metric is smooth away from the singularity, and one expects this regu- 
larity to persist for short time. Here we prove this. 

For r > 0, we let % r be the complement in § n+1 of the neighborhood of the 
North Pole in which ^j n i t < r. 

Lemma 11. For any n > and *i > 0, there exists C\ such that the Riemann 
tensor of any solution guj{t) is bounded by |Rm(x, *)| < C\ for all 

(x,t) e (S' i+1 x [*!,**]) u(3C ri x [0,*J). 

Furthermore, all covariant derivatives of the curvature are uniformly bounded: for 
any k < oo there is C'k such that |V fc Rm| < Ck on (S I1+1 x [*i,**]) U (X ri x [0,**]) 
for all co > 0. 

Proof. The curvature bound for * > *i was proved in Section 6.2. The derivative 
bounds for * > t\ then follow from Shi's global estimates by taking g{t\/2) as the 
initial metric [29]. 
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Now we estimate Rm in 0C ri for < t < t\ , for sufficiently small t\ . On dOC ri we 
have V'init = f\, and thus, for small enough tu > 0, ipuj(x, 0) = V'init = r i- By (2.5), 
we therefore have 

H<^<fr? for O^^j^L. 

We choose ii < r^/4(yi + n) and estimate |Rm| on d% ri using the v — function. 
For any r± > 0, the functions v^fr, t) corresponding to the metrics g^it) will be 
bounded away from zero on the interval hr\ < r < |ri and for < t < ti, if ti is 
small enough. Being solutions of the nondegenerate quasilinear parabolic equation 
v t = ^[v], all derivatives of the are uniformly bounded in the smaller region 
\f 1 /2r% < r < •y/3/2ri. Hence K, L and |Rm| are bounded on d% ri for < t < ti. 
Applying the maximum principle to equation (6.3), and possibly reducing £i one 
last time, we find that |Rm| is indeed bounded on % ri x [0, £i]. 

Observe that all derivatives V fc Rm are initially bounded, and also that they are 
bounded in the collar region, because they can be expressed there in terms of v u (r, t) 
and its derivatives. Thus the derivative estimates follow from a modification of Shi's 
local estimates, which allow one to take advantage of bounds on the curvature of 
the initial metric [10, Theorem 14.16]. □ 

6.4. Constructing the solution. Now we show that a subsequence of our regu- 
larized solutions converges to the solution described in our main theorem, i.e. to a 
smooth forward evolution from a singular initial metric go- We do this in two steps. 
We first demonstrate subsequential convergence to a limit solution that has all the 
properties we want, except possibly at the north pole NP. Then we find coordinates 
in which this limit is smooth at NP for all t > that it exists. Together, these 
results complete the proof of Theorem 1. 

Lemma 12. Any sequence LUi \ has a subsequence Ui k for which the solutions 
g Ufh (t) converge in C£ c on (§ n+1 \ {NP}) x [O.tJ. 

Proof. Let ajj \ be any given sequence. For arbitrary r\ > 0, we have just 
shown that the metrics g LOi (t) have uniformly bounded Riemann curvature tensors 
on % ri /2 x [0,t*]. Hence there exists a subsequence 0Ji k such that the solutions 
g Ui (•) converge in C°° on % ri x [0,i*]. For proofs, see Lemma 2.4 and Hamilton's 
accompanying arguments in [14], or else Lemma 3.11 and the subsequent discussion 
in [9, §3.2.1]. 

A diagonalization argument then provides a further subsequence which converges 
in on (S" +1 \ {NP}) x [0,tj. □ 

To complete the construction, we show in the next Lemma that the apparent 
singularity at the north pole is removable. 

Lemma 13. Let g*(t) be the limit of any convergent sequence g Ui (t). Then there is 
a homeomorphism $ : S n+1 — ► § n+1 for which the metrics g(t) = ^*[g*{t)j extend 
to a smooth metric on §" +1 for all t > that they exist. 
The homeomorphism is smooth except at the north pole. 

Because g*(t) satisfies Ricci flow, the modified family of metrics g(t) also satisfies 
Ricci flow on the punctured sphere. But since the modified metrics extend to 
smooth metrics on the whole sphere, they constitute a solution of Ricci flow that 
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is defined everywhere on S n+1 x [0,f*], except at the north pole NP at time t — 0. 
Its initial value is <&*[<7 ], so it is the solution we seek. 

Proof. Choose any time t\ € (0,f*), and let the homeomorphism $ be such that 

g(h) = (dx) 2 + tp(x, il) 2 5can- 

Then, since v(r, t±) is a smooth function, one finds that g(tx) extends to a smooth 
metric at the north pole. 

We write V for the Levi-Civita connection of g[t) at any given time t £ (0,t*), 
and V for the connection at the fixed time t\. 

We will show that the metrics g(t) extend smoothly across the north pole by 
finding uniform bounds for the derivatives \7 k g(t). These bounds imply that for 
any pair of smooth vector fields V, W defined near the north pole, the function 
g(t)(V, W) and its derivatives are uniformly bounded near the north pole, so that 
g(t)(V,W) extends smoothly. The metric g(t) therefore also extends smoothly. 

To estimate the derivatives V fc <?(£), we start with the identity V k g(t) = and 
then estimate the difference between the connections V and V. This difference is, 
as always, a tensor field, so that we may write 

d 

V = V + A, A = A\ h — <g> dx 1 ® dx k . 

The tensor A vanishes at time t = t\, and it evolves by 

(6.5) d t A) k = -g^VjRa, + V k R je ~ V e R ]k }. 

Our construction of g(t) was such that all covariant derivatives of the Riemann 
curvature of g(t) are bounded on any time interval e < t < t* . In particular, the 
Ricci tensor is uniformly bounded, and hence the metrics g(t) are equivalent, in the 
sense that 

(6-6) C^(t)~ 9{tl) ~ m ~ 

for some constants C\(t) (which become unbounded as t \ 0). Because of this, 
uniform bounds for some tensor with respect to one metric g(t) are equivalent to 
uniform bounds with respect to any other metric g(t'), as long as t,t' > 0. 

Since \7iRjk is bounded, (6.6) and (6.5) imply that d t A is uniformly bounded, 
so that we have 

(6.7) \A\ <C(t)|i-ii|, 

in which the constants C(t) again deteriorate as t \ 0. 

Because Vg(t) = (V — A)g(t) = —A * g(t) implies that Vg(t) is bounded, this 
leads to a Lipschitz estimate for the metric g(t). 

Next, we estimate by applying the time-independent connection V to both 
sides of (6.5). This tells us that 

d t VA = Vd t A = V(d t A) -A* A. 

We have just bounded A, so the second term is bounded, while the first term must 
also be bounded since it is given by (6.5). Thus we find that 

(6.8) \VA\ < C(t)\t - ti| and |VA| < C(t)\t - ti\. 

Inductively, one bounds all higher derivatives, W k A, and W k A. 

To bound W k F for any tensor F we write this derivative as V k F = (V — A) k F. 
Expanding this leads to a polynomial in F, VF, • • • , V k F and A, \7A, • • • , V k ~ 1 A. 
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The bounds for V ' J A which we have just established then show that bounds for 
F, VF, • • • , V fe F imply similar bounds for F, VF, • • ■ , V fe F. 

In the special case where F = g(t), we have the trivial bounds W k g(t) = 0, so 
that all derivatives V k g(t) are bounded. As stated before, this implies that the 
metrics g(t) extend smoothly to S n+1 . □ 

Appendix A. Rotationally symmetric neckpinches 

Here we recall relevant notation and results from [1, 2]. Remove the poles P± 
from S n+1 and identify §™ +1 \{-P±} with (—1,1) x § n . In [1], we considered Ricci 
flow solutions whose initial data are smooth SO(n+ l)-invariant metrics of the form 

(A.l) g = <p(x) 2 (dx) 2 + il>(x) 2 5can , 

where x G (—1,1). Parameterizing by arc length s(x) = J* (p(y) dy and abusing 
notation, we wrote (A.l) in a more geometrically natural form, 4 namely 

g= (ds) 2 + 4>(s) 2 g c 



/can ■ 



Smoothness at the poles is ensured by the boundary conditions that lim^^-ti ip s = 
=Fl and that ip/( s ± — s) be a smooth even function of s± — s, where s± = s(±l). 

Metrics of the form (A.l) have two distinguished sectional curvatures: let K 
denote the curvature of the n 2-planes perpendicular to a sphere {x} x S™, and let 
L denote the curvature of the Q) 2-planes tangential to {x} x S". These sectional 
curvatures are given by the formulas 

(A.2) K=-^ and L = ^f. 

Evolution of the metric (A.l) by Ricci flow is equivalent to the coupled system of 
equations 



(A.3) 



i> t = 4, ss - (n - l) 1 —^ = —tp[K + (n - l)L] 

W 

m/Jss „ 
tp t = — —(f = -nKtp 



in which d s is to thought of as an abbreviation of tp(x, t)~ l d x . 

In [1, 2], we called local minima of tjj necks and local maxima bumps. We called 
the region between a pole and its closest bump a polar cap. In [1], neckpinch 
singularity formation was established for an open set of initial data of the form 
(A.l) on §™ +1 (n > 2) satisfying the following assumptions: (1) The metric has 
at least one neck and is "sufficiently pinched" , i.e. the value of ip at the smallest 
neck is sufficiently small relative to its value at either adjacent bump. (2) The 
sectional curvature L of planes tangent to each sphere {x} x S" is positive. (3) 
The Ricci curvature Rc = nK ds 2 + [K + (n — l)L]tp 2 g can is positive on each polar 
cap. (4) The scalar curvature R — 2nK + n(n — 1)L is positive everywhere. In [2], 
precise asymptotics were derived under the additional hypothesis: (5) The metric 
is reflection symmetric, i.e. ip(s) — ip(—s), and the smallest neck is at x = 0. 

To describe the asymptotic profile derived for such data, let To < °o denote the 
singularity time, let \&(o", r) = ip(s, t)/ \/2(n — 1)(T — t) be the "blown- up" radius, 



4 The choice of s as a parameter has the effect of fixing a gauge, thereby making Ricci flow 
strictly parabolic. 
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where a = s/^/Tq — t is the rescaled distance to the neck and r = — log(To — t) is 
rescaled time. Then the main results of [2] can be summarized as follows. 

Theorem 7. For an open set of initial metrics satisfying assumptions (l)-(5) 
above, the solution (§ n+1 ,g t ) of Ricci flow becomes singular at a time Tq < oo 
depending continuously on g^. The diameter remains uniformly bounded for all 
t € [0, To]. The metric becomes singular only on the hypersurface {0} x S". The 
solution satisfies the following asymptotic profile. 
Inner region: on any interval \a\ < A, one has 

a 2 - 2 1 

^(cr, t) = 1 H h o(— ) uniformly as r — * oo. 

Intermediate region: on any interval A < \a\ < B^/t, one has 



0-2 

^((T, r) = sJ 1 + (1 + o(l)) — uniformly as r — > oo. 
Outer region: for < |s| <C 1, there exists a function h such that 



i>(s,t) = [l + h(\s\,T -t)] V " ' 



2 ^TogR' 

where h(a, b) —> as a + 6 \ 0. 

To flow forward from a neckpinch singularity, one must show that these asymp- 
totics can be differentiated. 



Lemma 14. Under the hypotheses of Theorem 7, the limit 
exists and satisfies 



ip(s) = lim ip(s(x,t),t) 
tST 



(A.4) = [1 + h (s)} V " ~ ' 



2 ^J— log s 

This asymptotic profile may be differentiated. In particular, one has 



(A.5) M*) = [l + hi(s)\ ^" ' ' 



2 ^J- log s 
for s > 0, where ho(s) — > and /ii(s) — > as s \ 0. 

Proof. Existence of the limit and equation (A.4) follow directly from the results in 
Section 2.19 of [2]. 

Recall that [2] proves there exist T» < T and x* > such that ip s (s(x,t),t) > 
and ip ss (s(x, t),t) > in the space-time region [0, x*] x [T*, T ). Moreover, the limit 
s(a;, To) exists for all a; e [0, x*). Fix any x, a; £ [0, cc*] with a; < x. Then one has 

~{s(x,i)-sfet)} = -n / K(s(x,t),t)ds =n ™>> ' j ' j ds > 

for all £ e [T,,T ), and 

fs(x,t) ps(x,T t ) 

ip(s(x ) t),t)-ip(s(x,t),t)= i[> a (s(x,t),t)ds > %/j s (s(x,t),t)ds >0. 

J s{x,t) J s(x,T,) 

Letting t /" T , one sees that ip(s) = ^>(s, T ) is monotone increasing in small s > 0, 
hence is diffcrcntiable almost everywhere. 
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In what follows, hi(s) denotes a family of functions with the property that 
hi(s) — > as s \ 0. First observe that 



(A.6) 



(-logs)" 1 / 2 ds 



V^logs 2 
[l + h 2 (s)\ 



V- log 1 ' 



Now suppose there exists a fixed e > such that 



Ma) > (1 + e) 



v^T 1 



2 v 7 - l°g s 

for a.e. small s > 0. Then applying (A. 4) and (A.6), one obtains 

^(S)= / s ^ («) rf g > (1 + s) [1 + fe 2 (g)] ^ ^— , 
Jo 2 V- log s 

> (1 + e) [l + ft 3 (a)]V(«). 

Dividing by ip(s) > yields 1 > (1 + e) [1 + h^s)], which is impossible for small s. 
It follows that 



M«) < [1 + /14(a)] 



Vn^l 1 



2 v^log^' 
The complementary inequality is proved similarly. 



□ 



Appendix B. Expanding Ricci solitons 

In this appendix, we prove a pair of lemmas which establish the assertion, made 
in Section 4, that every complete expanding soliton corresponding to a solution of 



(B.l) 



U 2 U PP + 

' P 



emerges from conical initial data, 



with [/„„ > 0. 



u 2 



Lemma 15. Let U(p) be a solution of (4.13) on an interval a < p < A < oo. 
// U' < and < U < 1, then lim p ^ A U{p) > 0. 

Proof. Since U is decreasing, we may assume that linip^^ U = in order to reach 
a contradiction. Assuming, as we may, that a is sufficiently large, we have 



IT 
U 



< 

= UU" 
< UU" 



\p- 



El 
p 

n-l 
n-l 



n — 1 



U 



a-u 2 ) 
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for all p > a. Integrate from a to b <G (a, A) to obtain 

,6 v , 



log U(a) - log U(b) 



U 



dp 



< 



n- 1 



UU 

= [uu'] h a - 

< U(b)U'(b) 

T) 

< -U(a)U'(a) + - 



dp 



(U') 2 dp + 

U{a)U'{a) 
1 



dp 
-1 



dp 



Keeping a fixed, we see that log U(b) remains bounded from above for all b £ (a, A), 
and hence that U(b) is bounded away from zero as b f A. □ 

Lemma 16. Let U(p) be a maximal solution of (4.13) defined for < p < A. 
Assume that < U(p) < 1 for all p G (0, A) and that 

lwaU(p) = l. 

Then A = oo, while U is strictly decreasing and is bounded from below by the 



positive quantity U a 



lim„ 



U{p). 



Proof. At any point where U' — 0, the differential equation (4.13) forces U" < 0. 
Thus every critical point of U is a local maximum. Since the solution stays between 
and 1 and starts at U(0) = 1, it can have no critical points for < p < A, hence 
must be decreasing. 

If the maximal solution exists, then cither A < oo and Yith p ^a U{p) = 0, or else 
A = oo. Lemma 15 rules out the possibility that U reaches at some finite value 
of p, so we must have A = oo. 

But when A = oo, Lemma 15 still applies and guarantees that [/(oo) > 0. □ 

Appendix C. The Bryant steady soliton 

A time-independent solution of the pde (4.17) satisfied by a solution V(a,t) in 
the inner region is a solution of the ODE £F[V] = 0. In this appendix, we find and 
describe all complete time-independent solutions. 

It is convenient here to consider U = VV . It follows from (2.7) that a stationary 
solution U must satisfy 



(C.l) 



d 

da 



dU 
da 



1 



(u- 1 u) 



ndU 
a da 



0. 



We introduce a new coordinate 

(C.2) ( = loga, 

in terms of which the metric can be written as 

(C.3) g - 



2 

d 3can 



9c 



a jj 2 scan - L jj2 

If one writes equation (C.l) in terms of this new coordinate, and defines 

TilJptn-lKr'-tO, 
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one sees that solutions of (C.l) correspond to solutions of the autonomous ODE 
system 

(C.4) d ^ = ( n +l)T + n(n-l)(U- 1 -U), % = T + (n - l)^" 1 - U). 

Observe that (T, U) = (0, 1) is a saddle point of this system. In fact, the lineariza- 
tion at (T, U) — (0, 1) is the system 

fV\ _ (n + 1 -2n(n - 1)\ ft\ 
\U>) ~ { 1 -2{n-l)) [uj 

with eigenvalues — (n — 1) and 2. The unstable manifold W u of this fixed point 
consists of two orbits. One orbit lies in the region U > 1,Y > 0. As £ - > oo 
this orbit becomes unbounded, and one finds that U — > oo. One finds that it 
does not generate a complete metric on K™ +1 . The other orbit lies in the region 
{(T,[/) : < U < 1,T < 0}. 




Figure 3. The (T, U) phase plane (n = 2) 



Lemma 17. The Bryant steady soliton is, up to scaling, the unique complete so- 
lution of (C.l) satisfying < U < 1 . 

Proof. It is well known that the Bryant soliton is the unique (up to scaling) com- 
plete, rotationally symmetric steady gradient soliton on for all n > 2. The 
point of the lemma is to exhibit it as a solution of (C.l), i.e. an unstable trajectory 
of (C.4) emerging from the fixed point (T, U) = (0, 1). 

We write the Bryant soliton (IR™ +1 , g, grad/) as it appears in [9, Chapter 1, 
Section 4], following unpublished work of Robert Bryant. The metric g is defined 
in polar coordinates on E" +1 \ {0} w (0, oo) x S" by 

(C.5) g = ds 2 + w 2 (s)g can . 

The soliton flows along the vector field grad/, where / is the soliton potential 
function, i.e. the solution of Rc(<?) + VV/ = 0. For this to hold, it is necessary and 
sufficient that w(s) and f(s) satisfy the system 

1 — w w 
(C.6) w ss = w s f s + (n - 1) f ss = n—. 
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Because U > for solutions of interest, we substitute a — w(s) and U — w s (s), 
thereby transforming (C.3) into (C.5). Then recalling from (C.2) that 

d( = — = — ds, 
a w 

one sees that (C.4) becomes 

T 1 - w 2 

(C.7) w ss =w s h(n— 1) T s = ra s / s + nro ss . 

w w 

The choice T = wf s transforms (C.7) into (C.6). 

In [9], following Bryant's work, the Bryant soliton is recovered from a careful 
analysis of trajectories of the ODE system [9, equation (1.48)] corresponding to 
x = U and y = nil — T near the saddle point (x,y) = (l,n) corresponding to 
(17, T) = (1,0). There it is shown that there exists a unique unstable trajectory 
(corresponding to U — * and T — ► — oo as s increases) that results in a (non- 
fiat) complete steady gradient soliton. Moreover, this solution is unique up to 
rescaling. □ 

Let U : (0, oo) — > (0, 1) be a solution of (C.l) which corresponds to the Bryant 
steady soliton. Every other such solution of (C.l) is given by U(ka) for some k > 0. 

We now note a few simple facts about the behavior of the Bryant steady soliton. 
Its properties near infinity are well known, but its properties near the origin are 
not as readily found in the literature. 

Lemma 18. 

(1) U is strictly monotone decreasing for all a > 0. 

(2) Near a — 0, U{a) 2 is a smooth function of a 2 , with the asymptotic expan- 
sion 

U(a) 2 = 1 + b 2 o- 2 + ^-b 2 ^ + g^L fria 6 + • • • , 
n + S (n + 3)(n + 5) 

where 62 < is arbitrary. 

(3) Near a — +00, U 2 has the asymptotic expansion 

,,, .9 9 4 — n 9 a (n — 4)(n — 7) o 

n — 1 (n — l) z 

where c 2 > is arbitrary 



Proof. (1) It is well known (see e.g. [9, Lemma 1.37]) that the sectional curvatures 
of the Bryant soliton are strictly positive for all a > and have the same positive 
limit at the origin. By (A. 2), the sectional curvature if of a plane perpendicular to 
the sphere {a} x S™ and the sectional curvature L of a plane tangent to the sphere 
{a} x §" are 

K = — and L = — , 

a g 

respectively. Hence U'(a) < and U(a) < 1 for all a > 0. 

(2) We showed in Lemma 17 that U = w' . Because w is an odd function that is 

smooth at zero [17], U is an even function that is smooth at zero. Thus U 2 has an 

asymptotic expansion near a — of the form 

U(a) 2 = 1 + b 2 a 2 + 6 4 CT 4 + b 6 a 6 + ■■■ . 
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By part (1), we must have 62 < 0. (Arbitrariness of 62 corresponds to the invariancc 
of U under the rescaling a t— > fcer.) The remaining coefficients are determined by 
the equation 3i n [[/ 2 ] = 0, which implies that 

= [2(n + 3)64 - 2n6 2 ] a 2 + 4 [(n + 5)6 6 - (n - 1)6364] cr 4 + ■ • ■ . 

In particular, 

64 = 6 2 > and ^6 = - — -6264 < 0. 

n + 3 n + 5 

The claimed expansion follows. 

(3) Let £ = cr -1 . Then {/ has an asymptotic expansion near £ = of the form 

[/(a) 2 = c + c 2 i 2 + c 4 £ 4 + c 6 £ 6 + • • • . 

It is well known (see e.g. [9, Remark 1.36]) that C~ l ^fs < w(s) < C^fs for s 
bounded away from 0. This forces Cq = 0. One finds that ^[J/ 2 ] = if and only if 

$ 2 U 2 Ul e - ^(£f/f ) 2 - (n - 1 - 3J7 2 )£[/f + 2(n - 1)(?7 2 - t/ 4 ) = 0. 

This allows C2 > to be arbitrary but forces 

4 - n , , - 4)(n - 7) , 

c 4 = -c 2 and c 6 = , — c 2 . 

n — 1 [n — \y 

The claimed expansion follows. □ 



In function 25 which we have used extensively in this paper is by definition 
(C.8) «B(<t) = U(ka) 2 

where U is as above in Lemma 18, and where the constant k is chosen so that 
(C.9) «8(<7) = \+o(a- 2 ), ((7->oo). 

Lemma 18 implies that such a choice of k exists, and that 55 is thus uniquely 
defined. 
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